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COMPARING TWO STATISTICAL ENSEMBLES OF 
QUADRANGULATIONS: A CONTINUED FRACTION APPROACH 

ERIC FUSY AND EMMANUEL GUITTER 


Abstract. We use a continued fraction approach to compare two statistical ensembles 
of quadrangulations with a boundary, both controlled by two parameters. In the first 
ensemble, the quadrangulations are bicolored and the parameters control their numbers of 
vertices of both colors. In the second ensemble, the parameters control instead the number 
of vertices which are local maxima for the distance to a given vertex, and the number 
of those which are not. Both ensembles may be described either by their (bivariate) 
generating functions at fixed boundary length or, after some standard slice decomposition, 
by their (bivariate) slice generating functions. We first show that the fixed boundary 
length generating functions are in fact equal for the two ensembles. We then show that 
the slice generating functions, although different for the two ensembles, simply correspond 
to two different ways of encoding the same quantity as a continued fraction. This property 
is used to obtain explicit expressions for the slice generating functions in a constructive 
way. 


1. Introduction 

The study of planar maps has given rise in the recent years to a lot of remarkable enu¬ 
meration results. A particularly fruitful approach consists in taking advantage of bijections 
between maps and tree-like objects called mobiles. This technique, initiated by Schaeffer 
[mi5] (reinterpreting a bijection by Cori and Vauquelin [7]) was extended in many different 
directions mm to deal with various refined map enumeration problems. Besides mobiles, 
another, slightly different, view on the problem consists in decomposing the maps into so- 
called slices, which are particular pieces of maps with nice combinatorial properties [5]. In 
particular, the generating functions for these slices were shown to obey discrete integrable 
systems of equations and in most cases, a solution of these equations could be obtained ex¬ 
plicitly. Moreover, the slice decomposition of a map is intimately linked to its geodesic paths 
and the knowledge of slice generating functions directly gives explicit answers to a num¬ 
ber of questions regarding the statistics of distances between random points within maps 

[a El m Ho]. 

A particularly important discovery was made in where it was shown that slice gener¬ 
ating functions happen to be simple coefficients in a suitable continued fraction expansion 
of standard map generating functions, making de facto a connection between the distance 
statistics within maps and some more global properties. On a computational point of view, 
this discovery provided a constructive way to obtain explicit solutions for the integrable 
systems at hand, by taking advantage of known results on continued fractions. 

Quite recently, the slice decomposition technique was used in m to describe the distance 
statistics of general families of bicolored maps, and, in particular, of bicolored quadrangula¬ 
tions, with some simultaneous control on the numbers of vertices of both colors. Explicit 
expressions for the corresponding bivariate slice generating functions were obtained in a con¬ 
structive way, leading in particular to explicit formulas for the distance dependent two-point 
function within bicolored quadrangulations. Remarkably, the expressions found for slice 
generating functions are very similar to those obtained (via a mobile formalism) in another 
problem of quadrangulations considered in [T]. There, the discrimination between vertices 
no longer relies on their color but rather on their status with respect to the graph distance 
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from a fixed origin vertex. Vertices namely come in two types: those, called local maxima 
which are further from the origin than all their neighbors, and the others. Bivariate slice 
generating functions can be defined so as to keep some independent control on the numbers 
of both types of vertices after the slice decomposition. Explicit expressions for these new 
bivariate slice generating functions were then guessed in [T] and, as just mentioned, their 
structure is very similar to that of their bicolored counterparts. 

The aim of this paper is twofold: first, we establish a strong connection between the 
problem of quadrangulations with a control on the vertex color, as discussed in [10) . and 
that of quadrangulations with a control on local maxima, as discussed in [T] . Then, we use a 
continued fraction formalism to re-derive, now in a constructive way, the explicit expressions 
found in [T]. 

The paper is organized as follows: in Sect|^ we introduce the two ensembles of quadran¬ 
gulations that we want to compare and define their generating functions at fixed boundary 
length. We then derive our first fundamental result which states that the fixed boundary 
length generating functions are in fact equal for the two ensembles. Sect. presents the 
slice decomposition of the quadrangulations at hand and shows that the corresponding slice 
generating functions may be obtained as coefficients of the same quantity, once expanded as 
a continued fraction in two different ways. In Sect.|^ we recall the integrable systems which 
determine the slice generating functions of both ensembles as well as the explicit solutions of 
these systems obtained in HO] and [T]. Sect. deals with results on continued fractions, and 
shows in particular how to extract their coefficients from those obtained via a standard series 
expansion. In one of the ensembles that we consider, the knowledge of this series expansion is 
not sufficient to get all the slice generating functions, a process which requires the knowledge 
of some additional quantity. An explicit expression for this latter quantity is conjectured in 
Sect. based on simplifications observed in the case of finite continued fractions, and we 
then show how it allows to recover the explicit formulas for the slice generating functions 
found in [T]. Sect. deals with another aspect of our problem, the existence of invariants, 
the so-called conserved quantities., as expected for discrete integrable systems. We show how 
to derive these invariants combinatorially for both ensembles and again emphasize the deep 
similarity existing between the conserved quantities for the two ensembles. We gather our 
concluding remarks in Sect. Some side results or technical derivations are presented in 
Appendices A and B. 

2. An equality between two bivariate generating fungtions for 

QUADRANGULATIONS WITH A BOUNDARY 

The aim of this section is to compare the generating functions of planar quadrangulations 
with a boundary weighted in two different ways, each of these weighting being bivariate, i.e. 
involving two independent parameters. As we shall see below, these two weightings, although 
fundamentally different, are intimately linked and some of the associated generating functions 
turn out to be equal. 

2.1. Two bivariate generating functions for quadrangulations with a boundary. 

Recall that a planar quadrangulation with a boundary denotes a connected graph embedded 
on the sphere which is rooted, i.e. has a marked oriented edge (the root-edge) and is such 
that all its inner faces, i.e. all the faces except that lying on the right of the root-edge, have 
degree 4. As for the external face, which is the face lying on the right of the root-edge, its 
degree is arbitrary (but necessarily even). As customary, the origin of the root-edge will be 
called the root-vertex. Let us now consider two particular different ways to assign weights 
to these maps. 

/First weighting: bicoloring the map. Since all their faces have even degree, planar quad¬ 
rangulations with a boundary may be naturally bicolored in black and white in a unique 


COMPARING ENSEMBLES OF QUADRANGULATIONS VIA CONTINUED FRACTIONS 


3 




Figure 1. Left: an example of rooted quadrangulation with a boundary of 
length 10, where each vertex is labelled by its distance to the root-vertex. 
The local maxima for this distance are indicated in gray. Right: the asso¬ 
ciated rooted general map with a bridgeless boundary of length 5, obtained 
by applying the Ambj0rn-Budd rule within each inner face (i.e. connecting 
the two corners within the face followed clockwise by a larger label) and, 
in the external face, connecting cyclically those corners followed by a larger 
label counterclockwise around the map. 


way, by assigning the black color to the root-vertex and demanding that no two adjacent 
vertices have the same color. We way then enumerate these quadrangulations by assigning 
a weight t, to each black vertex and a weight to to each white vertex. For convenience, the 
root-vertex receives a weight 1 instead of t,. We shall then denote by Fn = Fn{t,,to) the 
corresponding generating function for these maps with a boundary length 2n, i.e. with an 
external face of degree 2n. 

/ Second weighting: distinguishing local maxima of the distance. Our second weighting con¬ 
sists in giving a special role to the local maxima of the distance from the root-vertex. More 
precisely, we may label each vertex v of the quadrangulation by its graph distance d{v) from 
the root-vertex and look for the loeal maxima of this labeling, i.e. those vertices v having 
only neighbors with label d{v) — 1 (note that in all generality, neighbors of a vertex v may 
only be at distance d{v) — 1 or d{v) + 1 from the root-vertex). We decide to give a weight 
to to local maxima and a weight t, to the other vertices (see Fig. left). As before, the 
root-vertex receives a weight 1 instead of t, (note that the root-vertex can never be a local 
maximum). We shall call J„ = J„(t,,to) the generating function for these maps with a 
boundary of length 2n. 

The generating functions F„(t,,to) and J„(t,,to) may be understood as formal power 
series in t, and to , giving rise to convergent series for small enough , to • The first weighting 
is quite natural and was described in detail in [10] . We shall recall some of the corresponding 
results below. As for the second weighting, it may seem more artihcial but, as explained in 
[T], it arises naturally in two contexts: first, letting to —>■ 0 (i.e keeping the linear term in to) 
is a way to suppress local maxima of the distance, selecting quadrangulations arranged into 
layers between the root-vertex and a unique local maximum. These so-called Lorentzian or 
causal structures display a very different statistics from that of arbitrary quadrangulations 
[T]. As recalled below, the second weighting also arises naturally when enumerating general 
planar maps with a control on both their numbers of vertices and faces [Tj. 
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2.2. Equality of generating functions. Let us now prove a first fundamental equality, 
namely that 

( 1 ) Jn{t,,to) = Fn(t,,to) . 

To this end, let us recall the so-called Ambj0rn-Budd bijection of [1] between quadrangula- 
tions and general maps, slightly adapted to the case of quadrangulations with a boundary 
according to the rules of [J] . Starting with our quadrangulation with a boundary and label¬ 
ing each vertex v by its distance d(v) from the root-vertex, we associate to each inner face an 
edge as follows (see Fig.[^ right): looking at the corner^ clockwise around the face, exactly 
two corners are followed by a corner with larger label. We connect these two corners by an 
edge lying inside the original face. As for the external face, looking again at the corner labels 
clockwise around the face, i.e. counterclockwise around the rest of the map, exactly n corners, 
including the root-corner (lying immediately to the right of the the root-edge) are followed 
by a corner with larger label. We connect the n corners of this ensemble cyclically clockwise 
around the map, each edge connecting two successive corners in the ensemble (see Fig. [^. 
Finally, we mark and orient away from the root-vertex the edge connecting the root-corner 
to its successor. As explained in [1] (and its extension [4]), the obtained edges form a rooted 
planar map together with those vertices of the original quadrangulation which were not local 
maxima for the distance from the root-vertex. Each inner face of this map surrounds exactly 
one of the original local maxima, which get disconnected in the construction. More precisely, 
from mm, the above transformation provides a bijection between planar quadrangulations 
with a boundary of length 2n and rooted planar general (i.e. with faces of arbitrary degrees) 
maps with a bridgeless boundarj0of length n (i.e. with external face - lying on the right of 
the root-edge - of degree n and without bridge). The vertices of the quadrangulation which 
are not local maxima for the distance from the root are in one-to-one correspondence with 
the vertices of the general map while the vertices of the quadrangulation which are local 
maxima are in one-to-one correspondence with the inner faces of the general map. 

We may thus interpret Jn(t., to) as the generating function for rooted planar general maps 
with a bridgeless boundary of length n, weighted by t, per non-root-vertex and to per inner 
face. 

As for the label d{v) of a vertex v retained in this new map, it precisely corresponds to 
the oriented graph distance from the root-vertex to v on the new map, using paths oriented 
from the root-vertex to v which respect the following edge orientatiorQ all edges are oriented 
both ways except for the boundary-edges (i.e. the edges incident to the external face) which 
are oriented counterclockwise around the map. 

Forgetting about distances and labels, we may now use a standard construction to rebuild 
a quadrangulation with a boundary out of our general map. Coloring the vertices of the 
map in black, we simply add a white vertex within each inner face and connect it to all the 
corners within the fac^ By doing so, we get a bicolored quadrangulation with a boundary 
twice larger as that of the general map we started from, which we root by picking the edge 
leaving the root-vertex within the corner immediately to the left of the root-edge of the 
general map, and orienting it from its black to its white extremity (see Fig.right). Again 

'^Recall that a corner is an angular sector between two successive half-edges around a given vertex. The 
label of a corner is that of the incident vertex. 

^Strictly speaking, the extension |1] of [J shows that corners followed by a smaller label should be 
connected cyclically within all faces, including the inner faces, so that the resulting object is a hypermap, 
made of alternating black and white faces with, in our case, all black faces of degree 2 but one, of degree n, 
which we choose as external face. The Ambj0rn-Budd construction that we use here is recovered by squeezing 
all inner black faces, of degree 2, into simple edges while the external black face of degree n becomes the 
external face of the map. As for any face of a hypermap, its boundary is then necessarily without bridge. 

^In the underlying hypermap structure, the labels correspond to the distance using oriented paths going 
clockwise around the black faces. Squeezing the inner black faces of degree 2 results in simple edges oriented 
both ways, while the boundary-edges remain oriented oneway only. 

^Note that it is crucial that the boundary of the map be bridgeless for the obtained object to be connected. 
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Figure 2. Left: the rooted map on the right of Fig. with a bridgeless 
boundary of length 5. Right: The associated rooted bicolored quadrangu- 
lation with a boundary of length 10 obtained by inserting a white vertex at 
the center of each inner face of the map and connecting it to all the incident 
vertices around the face. 
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2n 


Figure 3. An example of directed path of length 2n = 12, made of elemen¬ 
tary steps with height difference ±1, starting and ending at height 0 and 
remaining (weakly) above height 0. The path is naturally colored in black 
and white. To compute we must sum over such bicolored paths with 
a weight Bi (resp. Wi) assigned to each descending step i —)■ i — 1 starting 
from a black height (resp. white height). 

this construction provides a bijection between rooted planar general maps with a bridgeless 
boundary of length n and planar quandrangulations with a boundary of length 2n, equipped 
with their (unique) bicoloration as defined in the previous section. The vertices of the general 
map are in one-to-one correspondence with the black vertices of the quadrangulation while 
the inner faces of the map are in one-to-one correspondence with the white vertices of the 
quadrangulation. 

We may thus interpret F„(t,,to) as the generating function for rooted planar general 
maps with a bridgeless boundary of length n, weighted by t, per non-root-vertex and to per 
face. Eq. 0 follows. 


3. Slice decomposition and continued fractions 


3.1. Slice decomposition for maps enumerated by As explained in [10], the quad- 
rangulations (with a boundary) enumerated by may be decomposed into slices by some 
appropriate cutting procedure. 

Labeling each boundary-vertex v by d(v), the sequence of corner labels, read counter¬ 
clockwise around the map starting from the root-corner, forms a directed path of length 2n, 
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apex apex 




Figure 4. A schematic picture of an i-slice contributing to Bi (left) and 
to Wi (right). 


made of elementary steps with height difference ±1, starting and ending at height 0 and 
remaining (weakly) above height 0 (see Fig. |^. Drawing, for each boundary-vertex v, its 
leftmost geodesic (shortest) path to the root-vertex and cutting along these geodesics results 
into a decomposition of the map into pieces, called slices. More precisely, to each descending 
step i —)■ f — 1 of the path corresponds an i-slice, defined as follows (see m for details): it 
is a rooted map whose boundary is made of three parts (see Fig. [^: (i) its base consisting 
of a single root-edge, (ii) a left boundary of length £ with 1 < £ < i connecting the origin of 
the root-edge to another vertex, the apex and which is a geodesic path within the slice, and 
(iii) a right boundary of length £—\ connecting the endpoint of the root-edge to the apex, 
and which is the unique geodesic path within the slice between these vertices. The left and 
right boundaries do not meet before reaching the apex (which by convention is considered 
as part of of the right boundary only). As a degenerate case when £ = 1, the left boundary 
may stick to the base, in which case the slice is reduced to a single root-edge. 

At this level, it is interesting to note that the distance dipf) from the root-vertex to any 
vertex v in the quadrangulation is directly related to its distance ds{v), within the f-slice it 
lies in, from the apex of this slice via 

(2) d(v) = ds{v) + i — £ 

if £ is the length of the left boundary of the slice. Indeed, it is clear by construction of the 
slices that either the root-vertex is the apex of the j—slice at hand or it does not belong to 
the slice at all and any path from v to this root-vertex must first reach one of the boundaries 
of the slice (possibly at the apex). In the first case, we have d{v) = ds{v) and * = so 
that (§ holds. In the second case, since the slice boundaries are part of geodesic paths to 
the root-vertex, d{v) is equal to ds(y) plus the distance from the apex of the slice to the 
root-vertex. In other words, d{v) — ds{v) has a constant value within the i-slice, which is 
obtained by taking for v the origin of the root-edge of the slice, namely d{v) — ds{v) = i — £, 
and ([^ follows. Note that, in an i-slice, i only acts as an upper bound on the length £ of the 
left boundary. The vertices v of an i-slice may then be labelled by non-negative integers in 
two natural ways: either by their distance ds{v) to the apex or by this distance plus i — £. 

Let us call Bi = Bi{t,,to) (resp. Wt = Wi{t,,to)) the generating function for i-slices 
whose root-vertex is black (resp. white), with a weight i, per black vertex and to per white 
vertex except for the vertices of the right boundary (including the endpoint of the root-edge 
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and the apex) which receive a weight 1 The slice decomposition implies that m 
(3) F„ = Zo+o(2n; , 

where ZQg(2n; denotes the generating function of paths of length 2n, made 

of elementary steps with height difference ±1, colored alternatively in black and white, 
starting and ending at black height 0 and remaining (weakly) above height 0, with each 
descending step from a black height i to a white height i — 1 weighted by Bi and each 
descending step from a white height i to a black height i — 1 weighted by Wt (see Fig. [^. 

The set of identities (|^ for all n > 0 can be summarized into the continued fraction 
expression 


(4) F(^)^^F„z" = 


n>0 


1-z- 


Wi 


1 - z- 


Bo 


1 - 


W, 


1 - Z: 


B4 


1 - 


with the convention that Fq = 0 and where F(z) = F(z;tm,to) implicitly depends on t, and 

to- 

3.2. Slice decomposition for maps enumerated by J„. Let us now play the same 
game with maps enumerated by J„, which are the same maps as those enumerated by F^, 
but now with the second weighting. We may again apply the same slice decomposition, 
resulting in the same i-slices as before. More precisely, labeling each boundary-vertex v by 
d{v) gives rise to a path of length 2n (from height 0 to height 0, remaining above height 0) 
and each descending step i —)■ i — 1 gives rise to an i-slice. To assign the second weighting 
to the quadrangulation, we must label each vertex v of the i-slice by its distance d{v) (in 
the quadrangulation) from the root-vertex of the quadrangulation. As explained above, if 
the 7-slice has a left boundary length £ {1 < £ < i), this amounts to label v by ds{v) + i — £ 
where dg {v) is its distance (within the slice) from the apex of the slice. To recover the correct 
weights, we must first give weight 1 to all the vertices of the right boundary (including the 
endpoint of their root-edge and the apex) in order to avoid double weightings after regluing 
the slices. As for the vertices lying on the left boundary of the slice and different from the 
root-vertex of the slice, they cannot, as part of a geodesic of the original quadrangulation, 
be local maxima as they have a neighbor with larger label along the geodesic path. They 
receive a weight t, accordingly. Considering now vertices v lying strictly within the slice, 
they have all their original neighbors lying in the slice and, from § , are local maxima for 
the distance d{v) if and only if they are local maxima for the distance ds{v) within the slice. 
For all these vertices, we may thus use the distance ds{v) within the slice to detect the local 
maxima, and give them the weight to, while non-local maxima for ds{v) get the weight t,. 
The last vertex to consider is the origin w of the root-edge of the slice: for this vertex to 
be a local maximum of the original distance d(v), it must both be a local maximum of the 
distance ds{v) within the slice and have no neighbor with larger label after regluing. Now two 
situations may occur: either the boundary-vertex preceding w along the boundary (oriented 
counterclockwise around the quadrangulation) has label d{w) — 1 and then the slice at hand 
sticks to the boundary so that w has all its neighbors within the slice. Then w is a local 
maximum for d(y) if and only if it is a local maximum for ds{v). Or the boundary-vertex 
preceding w has label (i(w)-l-l is which case w is not a local maximum for d{v), irrespectively 
of whether or not it is one for ds{v). 


®The fact that these vertices receive a weight 1 is to avoid double weighting upon re-gluing the slices into 
a quadrangulation. Indeed, all these vertices are already part of a left boundary, except for the root-vertex. 
In the end, only the root-vertex gets a weight 1, as wanted. 
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Figure 5. A schematic picture of an i-slice contributing to Qi (left) and 
to Pi (right). Local maxima of the distance from the apex are indicated 
in gray and non local maxima in black. The root vertex may be a local 
maximum or not, and receives the weight to or t, accordingly in Qi, while 
it always gets the weight in Pi. 



Figure 6. An example of directed path of length 2n = 12, made of elemen¬ 
tary steps with height difference ±1, starting and ending at height 0 and 
remaining (weakly) above height 0. To compute Jn, we must sum over such 
paths with a weight Qi (resp. Pi) assigned to each descending step i —)■ i — 1 
following an ascent (resp. following a descent). 

To summarize, we are led to consider two different generating functions for i-slices. In 
the first generating function Qi = Qi(t,,to), all the vertices of the f-slice receive a weight Q 
or t, according to whether or not they are a local maximum for the distance ds{v) from the 
apex within the slice (in particular the vertices of the left boundary different from the root- 
vertex of the slice get a weight t, as wanted), except for the vertices of the right boundary 
(including the endpoint of their root-edge and the apex) which receive a weight 1. In the 
second generating function Pi = Pilt,,to), we assign exactly the same weights as in Qi, 
except for the root-vertex of the slice, which gets the weight t, irrespectively of whether or 
not it is a local maximum for ds{v) (see Fig. [^. 

Returning to the slice decomposition, it now implies, for any positive n, that 

Jn = ^(yo(2”-; {Pi}i>lAQi}i>l) ) 

where ZQg(2n; {Pi\i>i,{Qi}i>i denotes the generating function of paths of length 2n, made 
of elementary steps of height difference ±1, starting and ending at height 0 and remaining 
above height 0, with each descending step from height i to height i — 1 weighted by Pi if 
it follows a descending step i + 1 ^ i and by Qi if it follows an ascending step z — 1 —>■ f 
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Figure 7. An example of directed path of length 2n = 12, made of elemen¬ 
tary steps with height difference ±1 and elongated steps (of horizontal length 
2) with height difference 0, starting and ending at height 0 and remaining 
(weakly) above height 0. To compute J„, we must sum over such paths 
with a weight ^ 21-1 (resp. ^ 21 ) assigned to each elongated step i — 1 —>■ f — 1 
(resp. each descending step z —>■ i — 1). 


(see Fig. [^. Setting Jq = 1 and using the shorthand notation J{z) = this is 

summarized into the new continued fraction expansion 


( 5 ) 


J{z) = Y, JnZ^ 

n>0 


1 


l-z(Qi-Pi) 


Pi 

1-z{Q2-P2)-z - — - p— 

^-ziQ3-P3)-Z 


which we may write as 


( 6 ) 


J(z) = 


1-zYi- z- 


Y2 


I-ZY 3 - z- 


Y4 


1 — zY^ — z 


Ffi 


1 - 


upon dehning 

(7) Y2,-i = Q,-P,, Y 2 , = P. 

for z > 1. To understand ([^, or equivalently Q, we note that, expanding the right hand side 
of this latter equation, the term of order z" builds the generating function ZQg(2n; {Fi}i>i) 
of paths of length 2n starting and ending at height 0 and remaining above height 0, made of 
elementary (i.e. of horizontal length 1) steps of height difference ±1 together with “elongated 
steps” of horizontal length 2 and height difference 0. Each elementary descending step from 
height z to height z — 1 (z > 1) receives a weight l 2 i while each elongated step at height 
z — 1 (z > 1) receives the weight ^ 21-1 (see Fig. [^. Deforming each elongated step at height 
z — 1 into a sequence of elementary steps z—1—>z—>z — 1, we recover paths made only 
of elementary steps of height difference ±1, and (after regrouping all paths with the same 
deformation) receiving a weight 121-1 + T 2 i = Qi for each sequence z — 1—^z—^z — lor 
equivalently for each descending step z —>■ z — 1 following an ascent, and a weight l 2 i = Pi 
for those elementary steps z —>■ z — 1 which are not part of a sequence z—1—^-z—^-z — 1, 
i.e. follow a descent. In other words, Zf^Q(2n; {Yi}i>i) = .Z^g(2n; {Pi}i>i, {(5i}i>i), which 
explains the identity ([^. 
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Figure 8 . A schematic picture of the slice decomposition of a (non-trivial) 
i-slice into two slices, leading to the recursion relations Q and (H). 


To conclude this section, let us rewrite our fundamental equality Q in the more compact 
form 

J{z;t,,to) = F{z;U,to) ■ 

4. Getting the slice generating functions by solving recursion relations 

The slice generating functions Bi, Wi, Pi and Qi satisfy systems of non linear recursive 
equations which may be derived by performing a slice decomposition of the slices themselves. 
Indeed, when the Tslice, of left boundary length i, is not reduced to a single edge, we may 
look at the sequence of vertices encountered clockwise around the face lying on the left of 
the root-edge of the slice and draw the leftmost geodesic paths from these vertices to the 
apex. Using the labeling ds{v) + i — £ , the sequence of encountered labels, starting from its 
root-vertex, is either i—>- 1 - 1 - 1 —>- 1 —>■* — l,i—>■*—1 —>■*—>- 1 —fori—>•*—1—— 2—1 
(if i > 2) and, upon cutting along the leftmost geodesic paths, a new slice arises for each 
descending step of this sequence (see Fig.[^. 

For the first weighting, this decomposition, applied to f-slices enumerated by Bi and Wi, 
leads to the system 


Bi — it + Bi(Wi-i + Bi -h Wi+i) 

W, = to + IU,(U,_1 + VF, + U,+i) 


for i > 1, with Bq = Wq = 0. For the second weighting, this decomposition, applied to 
i-slices enumerated by Pi and Qu leads similarly to the system 


Pi — t, + Pi{Pi-i + Qi + Qi+i) 

Qi = to + Qi{Pi-i + Qi) + PiQi+i 

for i> 1, with Pq = Qq = 0 . 

The solution of Q was derived in [10] . Parametrizing and to by x and 7 via 


U 


a ;(7 — a:)^(l — 


(10) 



2 


to 


x{'y — ai^)(l — 72:)^ 
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with Ixl < l[^ it was shown that 

^ ^^_(l-x2*)(l-7a:2*+3) 

2l -D- 

( 11 ) 


(1 — 7 a: 2 *+^)(l — 


W2^+l = W 


B2i+1 — B 
for i > 0 , where 
B = 


(1-^2»+1/7)(1-^2z+4) 

(1 — a: 2 *+ 2 )(l — a:^*+3/7) 
x(7 — xY 


Wo,. = W 


(1 - 'yx^^+^)(l - x^^+-^) 
(1 — X^^+^)(l — 7x2®+3) 
(1 — a;^®)(l — a;^®+^/ 7 ) 


(1 — a;2®+V7)(l — a; 2 *+ 2 ) 


W = 


X + x^ + j — 6 x^j + x'^j + xj^ + x^j^ 
x(l — ^xY 


a: + + 7 — 6 a ;^7 + a :‘^7 + x 7 ^ + ^^ 7 ^ 

As for the solution of ([^, it was guessed in [T]. Parametrizing now t, and to by y and a 


via 


( 12 ) 


t, = 

to = 


2/(1 - ayYY - oty^) 


(1 + y + ay — 6 ay^ + ay^ + a^y^ + a^y"^)^ 

a2/(l - yfiX - oi^yY 


(1 + y + ay — 6 ay^ + ay^ + a^y^ + a^y^)^ ’ 
with |y| < 1 , it was found that 


(13) 


P,=P 


for i > 0 , where 
P = 


(1 - 2/')(l - Qy"+ ) 

(l_y*+l)(l_Q,yi+2) 

2/(1 - ay? 


Q, = Q 


(1 - y")(i - a^y"+^) 

(1 — ay*+^)(l — ay*+^) ’ 


(14) 


Q = 


1 + y + ay — 6 ay^ + ay^ + a^y^ + a^y'^ 

_ay(l - yf _ 

1 + y + ay — 6 ay^ + ay^ + a^y^ + a^y^ 


Note that the two parametrizations (10 1 and (12) are actually equivalent providing we relate 
y and a to x and 7 via 


1 


7 


2 ’ 


y = 7 x 


With this correspondence, we immediately deduce that 
P = B , Q = W. 


This should not come as a surprise since, from (11) and (13), B, W, P and Q are the i —>■ 00 
limits of Bi, Wi, Pi and Qi, enumerating slices with no bound on their boundary lengths. 
From Q and Q, both pairs {B, W), and (P, Q) are determined by the same closed system^ 
namely: 

B = U+B{B+ 2W) , IF = to + VF(1F + 2P) , 

P = t, -\- P{P + 2Q) , Q = to Q{Q + 2P) . 

Let us end this section by rewriting the results for Pi and Qi in terms of F, as defined in 
0. First, Eq. § may be rewritten as 

Ui = t, + F2i(F2i-2 + F2i_i + Y2i + Y2i+1 + F2i+2) 

F 2 i_l = {to — t.) + F2i-l(F2i_2 + F 2 i_i + Y 2 i) 


for i>l, with Ig = 0. From (13), we immediately deduce the solution 

(1 - 2 /*^^)(l - ay*+^) 


(17) 


Y2,=P 


(1 - 2/')(l - ay"~^Y 

(l_y*+l)(l_Q,yi+2) 


F 2 *+i = F 


(l_y*+ 2 )(l_ayi+ 2 ) ’ 


®The parametrization in invariant under ( 31 , 7 ) I/ 7 ) so we may always choose |fc| < 1. 
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for i > 0, with 

Y=Q-P. 

The aim of this paper is to go beyond the guessing approach of [T] and to provide a con¬ 
structive way to obtain this latter formul a (|17[ ), and consequently 0, upon using general 
results for continued fractions of the type (1^. This is indeed the constructive approach used 
in [To] to obtain the expressions 0 from general results for continued fractions of the type 


5. Getting the slice generating functions by extracting continued fraction 

COEFFICIENTS: GENERALITIES 

5.1. The Stieltjes type. Eq. Q is a continued fraction of the so-called Stieltjes type. Its 
coefficients B 2 i and W 2 i-i for * > 1 are known to be related to the coefficients via the 
relations 


(0) ,;,(!) 


(>«) 


Li'j; / l(i) 
'H-l 'b-2 


W2^-l = 


''■i—2 'h—2 


for i > 1, in terms of the Hankel determinants 

— det(Tl2_|_^)o<n,m<2 = det(Tl^_|_j72+l)o<n,m<i 

for z > 0, with the convention = 1. These expressions were used in [TU] to 

obtain the expressions ( |TT| for B 2 i and W 2 i+i (i > 0). As for the expressions of B 2 i+i 
and W 2 i, it is clear from (j^ that B^ and play symmetric roles upon exchanging t, and 
to ■ The expressions 0 for i? 2 i+i and W 2 i are simply deduced upon this transformation, 
which amounts to a change yol/y, BoITin the formulas (see [IHI)- At this stage, it 
is important to note that the knowledge of the generating functions Fn is not sufficient to 
determine all the BiS and Wi’s as the associated continued fraction involves only one parity 
of the index i {Bi’s with even index i and WiS with odd index i) and that we have to rely on 
a symmetry principle to get the other parity. Otherwise stated, the derivation of all the BiS 
and WiS requires in principle the knowledge of a second family of generating functions. In 
the present case, these generating functions are nothing but those of rooted quadrangulations 
with a boundary of length n, bicolored in such a way that their root-vertex is white instead 
of black. Of course, by symmetry, those are nothing but the n > 1 and a simple 

symmetry argument is sufficient to conclude. 


5.2. The type of Eq. ([^. When dealing with a continued fraction of the type of Eq. §, 
a first remark should be emphasized: the knowledge of J„ is not sufficient to determine the 
coefficients Yi. Indeed, expanding in z gives rise to the first equations: 


Ji — (Yi -I- T 2 ) 

(19) -^2 = (W + Y 2 f + Y 2 iY 3 + 1^4) 


and it is easily seen that, at each step, two new T^’s appear on the right hand side, so that 
the system is clearly underdetermined. 

As shown in lillj, a full determination of the coefficients Yi requires, in addition to the 
set of J„ for n > 1, the knowledge of IT and of a second family of quantities J„ = J„(t., G), 
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n > 0, satisfying 
(20) J{z) = Y, JnZ^ 

n>0 


1 


l-zYi- z- 


r2 


l-zY.- z- 


Y. 


1 — zY^ — z- 


Y. 


1 - 


where we have defined (assuming h) 7^ 0 for alH > 1) 


(21) f2z-l = 


1 


Y2i = 


Yo, 


\r \r 

^2i-l -i22-1^2i+l 

for i> 1. Expanding in z now gives rise to the first equations: 


( 22 ) 


7 (1^2 + Fa) 

j _ {Y2 + Y^f 
' (ldi3)2 


1^2 (^4 + ds) 
YiYiY^ 


Knowing Yi, the first equation of (19) yields Y 2 , then the first equation of (22) yields Ya, 


the second equation of (19) yields Yj, and so on. The Yj’s are now fully determined and a 


compact formula may be written as follows: define 
1 if n = 0 


(23) 


Jn 


= Yi Jn—l 


J-r 


if n > 1 


if n < —1 


and the Hankel-type determinants 
(24) 


= det{j„+m-i) 


l<n,m<i 


Then we have, for z > 1, the following formulas, reminiscent of (18), 


(25) 


Y2^ = 


H, 


(0) 

i-l 




( 1 ) 


H 


(0) 


H, 


( 1 ) 

i+l 


^2.-1 = ^ 


■h: 


(0) 


H. 


(0) 

i-l 


with the convention = Hq 


( 1 ) 


= 0. A proof of these formulas can be found in PE]. 
We present a slightly simpler proof in the Appendix A below. To summarize, when dealing 
with a continued fraction of the type of Eq. (§, we may extract the coefficients Y^ if, in 
addition to J(z), we also know Yi and J{z). As we shall see in Sect. below, getting a 
simple expression for Yi combinatorially may be achieved upon using a so-called conserved 
quantity. As for J(z), we have not been able to obtain it via combinatorial arguments 
(as opposed to the previous section, we cannot rely here on any symmetry principle to get 
from Without the knowledge of J{z), Eq. ([^ yields a much weaker 

system than the recursion equations (16). In fact, any arbitrary choice of Jn will lead, 


through (25), to a set ofYi’s satisfying Eq. (|^, while the actual Y^’s, solution of Eqs. (16), 
correspond to a unique value of the J„’s, to be determined. 

As we shall now explain, we may however conjecture a simple expression for J(z), based 
on an explicit solution of the problem in the case of finite continued fractions. With this 
conjectured form of J(z), we may then verify that the obtained Y^’s precisely match their 
actual expressions ^ guessed in [T]. 
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5.3. The case of finite continued fraction. In this section, let us briefly digress from our 
combinatorial problem and discuss the case of a finite continued fraction. More precisely, let 


J{z) 


1 


1-zYi- z - 

I-ZY3-Z 


Y 2 


Y 4 


1 - zY2a-3 - Z- 


Y2a-2 


^ — Z Y2a-1 

where hi, I2, ■ • • , Y 2 a-i denote independent indeterminates. We also define 
J{z) = - ^ 


1-zYi- z- 


Yn 


1-zY.- z- 


14 


1 - zY2a-3 - Z- 


Y2a-2 


1 - zY2a-l 


and Yoi = 


To, 


for 1 < z < a . 


with T2i-i =47^ for 1 < z < a - 

r2i-l 421-1421+1 

The rational function J{z) is easily seen to be the ratio of a polynomial of degree a — 1 in z 
by a polynomial of degree a in z, hence characterized by (a—1 + 1) +(a+ 1) — 1 — 1 = 2a — 1 
coefficients (the last two — I’s correspond to removing a global factor in both the numerator 
and the denominator, and ensuring that J(0) = 1), depending on the 2a — 1 indeterminates 
Ti, +2, • • • , T2 q_i. In this case, the knowledge of the function J(z) alone therefore entirely 
determines all the coefficients of the continued fraction. This property may be reconciled with 
the apparently contrary statement of the previous section by noting that, in the present case 


of a finite continued fraction, both Ti and J(z) (defined via (20) and (21)) can be deduced 


from J(z). More precisely, we have the following relations, derived in Appendix A below: 

1 


(26) j{z) = - — j(- 
z \ z 


Ti = - 


lim zJ{z) 

Z—¥(X> 


Note that J{z) is also a rational function of z and that the expression for Ti simply rephrases 
the desired property that J(z) = 1 + 0{z). Knowing J(z), Ti and J{z), we can then deduce 
the coefficients for all integer n via their definition (23) and get Y 2 ,Y^, ■ ■ ■ ,Y 2 a-i from 


Eqs. (24) and (25) (which are also valid in the case of a finite continued fraction - see 


Appendix A). The relations (26) are proved in the Appendix A below. 


6. Recovering (13) from the continued fraction formalism 


6.1. A conjectured expression for Ti and J{z). Returning now to our enumeration 
problem of z-slices with the second weighting, let us conjecture that, although our continued 
fraction is now infinite, the relations (26) still hold for the particular choice of Yi we are 
interested in, namely the solution of (16). More precisely, J{z), originally defined as a 
power series in z, is convergent for small enough real z (namely 0 < z < l/(-\/Q + VP)^, 
see explicit expressions below - here we assume that t, and to are small enough positive 
reals so that P and Q are positive reals) but may be analytically continued to large enough 
real z (z > l/{y/Q — ^/Pf^). This allows us to define Jiftjz) for small real z (namely 
0 < z < (>/Q — '/P'f') s-iid our conjecture is that, in this range, J(z) is obtained via the 
relation J(z) = —(Ti/z) J(l/z) with a value of Ti adjusted so that J(0) = 1. Assuming this 
property, let us now see if we can then recover the desired expression (13), or equivalently 


0- 
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Figure 9. A sketch of the three interpretations of the generating function Z(^Q(2n; P, Q). 


Let us start by recalling the expression of F{z), hence J{z). From [TU], we know that 
(27) F„ = ^(1 - B - VF) Z+o{2n; Z+^{2n + 2; B, IF) 


t. 


where ZQQ{2n-,B,W) denotes the generating function of paths of length 2n, made of ele¬ 
mentary steps with height difference ±1, colored alternatively in black and white, starting 
and ending at black height 0 and remaining (weakly) above height 0, with each descending 
step from a black height to a white height weighted by B and each descending step from a 
white height to a black height weighted by W. A derivation of this expression via slices is 
recalled in Sect. 0 below. 

Equivalently, since Jn = Fn, P = B and Q = W, we have 
Jn = Aq Z^Q(2n; P, Q) + Ai ZQ Q{2n + 2; P, Q) , 

(28) p p 

Ao = -(l-P-Q), A^ = --. 

Let us introduce 

Z(z;P,g) = ^Zo+o(2n;P,Q)z" , 

n>0 

which, by definition, is a solution of 

1 


(29) Z{z;P,Q) = 


1-z- 


Q 


l-zPZ(z;P,Q) 

Note that this (quadratic) equation in Z is equivalent to the equation 

1 


Z{z-P,Q) = 


l-z{Q-P)-zPZ{z-P,Q) 
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so that Z(^Q(2n; P, Q) is also the generating function of paths of length 2n, made of elementary 
steps with height difference ±1, with each descending step weighted by Q if it follows an 
ascending step and by P otherwise. Alternatively, Z^g(2n; P, Q) enumerates paths of length 
2 n, made of elementary steps of horizontal length 1 and height difference ±1, and elongated 
steps of horizontal length 2 and height difference 0, each elongated step receiving the weight 
Y = {Q — P) and each elementary descending step the weight P (see Fig. [^. As a continued 
fraction, we thus have 


(30) Ziz;P,Q) = 


1 


l-zP-z- 


P 


Y = Q-P. 


l-zY-z~ 


P 


1- zY - z- 


1 _ ., 


Z{z-,P,Q)-l 


In terms of Z, we may write 

J(z) = Aq Z(z\ P, Q) + Ai 

z 

and, in components 

(31) = Ao [z^]Ziz- P, Q) + Ai [z"+i]Z(z; P, Q) 


for n > 0. At this stage, it is important to note that Eq. (29) yields two branches for Z for 
real z, namely 

1-Y z- ^{l-Y zf -APz 


Z_{z-P,Q) = 
Z+{z-P,Q) = 


2Pz 

l-Yz+^{l-Yz)'^-APz 

2Pz 


Y = Q-P 


for 1^1 < l/{\/Q + VP)^- To recover the coefficients J„, we must expand J{z), hence 
Z{z] P, Q) at small z, which requires to choose 

Z{z;P,Q) = Z_{z-P,Q). 

From Z_{z] P, Q) = l-\-Qz + 0{z‘^) we get J{z) = l + 0(z) (since from (15), Aq-\-AiQ — 1), 
as wanted. Using ( [^ , we find that 

J{z) = -Fi (Ao Z{z- P, Q) + Ai (z Z(z; P, Q) - 1)) 


1 


1 


Z+(z;P,Q) = 


Y = Q-P 


where Z[z] P,Q) = -Z P,Q 
Again, we have two possible branches for real z: 

■-F + V(z-F)2 -4Pz 
2 Pz 

and, to get the J„’s, we must, depending on whether P > Q or Q > P, choose the first or 
second branch respectively to get rid of the 1/z term when z —>■ 0. Both situations yield 
actually the same expression for J„. Assuming for instance P > Q, we get 

J„ = -yiAo [z"]Z_(z;P,Q)-yiAi [z"](zZ_(z; P, Q) - 1) . 

The value of Yi is obtained by ensuring that Jq = 1. Using F_(z; P, Q) = —l/Y + 0(z), we 
deduce 


(32) ( ^ + Ai ) = 1 , hence Yi = 


Y 


Aq + Ai Y 


(Q-P)(1-P-2Q) 

i-2g 


This value matches that obtained directly via conserved quantities in Sect.[^ 
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Using the identity 

Z-iz;P,Q) = j^z. 


,2-^p^Q 


we eventually arrive at 
Yi 


Jj). — 


(33) 


(P-Q)2"+i 
Yi 


{P-QY 

{Ao [z^]Z_iz;P,Q) + A, (P - Q)"[z"](zZ_(z; P, Q) - 1)) 


{Ao Z+,(2n; P, Q) + A^iQ - Pf Z+,(2n - 2; P, Q)) 


(Q - P)2"+i 

for n > 0. For Q > P, we must use instead Z^{z\ P, Q) but again in this case, ^+(z; P, Q) = 


-l/Y + 0{z) and Z+{z\P,Q) = pYgZ^ 


(p_Q)2,-r, 


pq) 


so that the expressions (32) and 


(33) remain unchanged. The first line of ( |33[ ) may be rewritten as 

(Ao [z^]Z. (^;P,q) + [2-] (zZ. (^;P,q) - l)) 

(Ao lz"]Z(z;P,Q)+A^ [z"-i]%;P,Q)) 


(34) 


J - ^ 

'Jn — y 

Y 


for n > 1, with 

Z(z; P, Q) = 




1 - zF- 


P 


1-zY-z- 


P 


l-zY-z 


P 


1 - 


1 ~ P 

where Y = — and P = . 

Y y 2 


Upon defining 


kn = 


1 if n = 0 

Y Zn-i if n > 1 


Z-n if n < —1 

where = [z'^\Z{z-, P,Q) = Zo+o(2n;P,Q) and = [z"]F(z;P,Q) = Z+^{2n-P,Q) jY^'^ 
(n > 0), we may summarize (31), (32) and (|34| into 


Y Y 

jn = Ao Y f^n + Ai y kn +1 for all integer n. 

6.2. Computation of and hY\ The above expression for j„ for all integer n opens 
the way to compute and via (24). Indeed, as shown in Appendix B, the coefhcients 
kn satisfy a set of linear relations of the form 


X', 


(i-i) 


b-i) 


= 0 , 


(35) ^ kn-^+mi-iy 

m=0 1 

while, for n = 1 and n = i + 1, we have 


2 < n < i 


Y, ki.,+^i-iy 


(i-i) 


31 — 1 


( 36 ) 


m—0 


X 


(2-1) ^2(2-1) 


2-1 


E ki+m{-iy 


P 

I i—l — m ^J 


m—0 


X 


( 2 - 1 ) 

2-1 


Yi 
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From these relations, replacing the first column {jn-i)i<n<i of by the linear combination 

i-l ~ ~ 


m—O 

us to write 


of this first column with the i — 1 last ones allows 


l<n<i 


( 0 ) = 


^0 


Fi P 


i-l 


Y Y2(i-1) 
0 






Y y* 


— 2-1-2 

Jo 

-2-1-3 

Jl 

Jo 

ji-2 

Jl 

ji—l 






Alternatively, the coefficients fc„ satisfy another set of linear relations of the form (see Ap¬ 
pendix B) 

i-l 

(37) ^ =0 , 2<n<i 

m—0 

while, for n = 1 and n = i + 1, we have 

-Di—l 


^ = (- 1 ) 


(38) 


m—0 


pt 

i-l ^ 


^ fc.+i_^(-i)™4*-i) = yp*-i(y + p). 


m—0 


From these relations, replacing the last column {jn)i<n<i of by the linear combina¬ 


tion E in-m(-l)^ 


m—0 


of this last column with the i — \ first ones, and using 


l<n<i 


Xg* = 1 (see Appendix B), allows us to write 




J-i+2 

j-i+3. 

Jo 

Jl 


JO (-ir^^oy^ 

Jl 0 


To summarize, and are fully determined by the system 

31—1 pi—1 

^ Yi 

ji-lt 


Ji-2 


J-_1 AiFrp*-i(y+p) 
r(i) 


p 


>i— 1 


= ^ 0^1 y^Hll\+A^Y,P^-\Y+P)Hll\ . 


r 7 "( 0 ) _ A V ^ I A W ^ 

^ + A,Y,P^-^{Y + P)Hf^\ 
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i •! 

) 

2 — 1 1 

[ 

i 

1 



i-2> 

1 

cc 

o * 


CO 

0) < 


'a _ 

’P - 


hard piec 

+ 

2 i 

2 ' 


2 ' 

1 i' 

1 ' 

, 1 

1 < 


Figure 10. A schematic picture of Eq. (411, identifying as the gener¬ 
ating function of hard pieces on a linear graph with i — 1 vertices. 


for i>\ with ' = Hq = 1. Upon setting 

R»-i) 

/V\ 

(39) 


l‘»> . I ^ 


H. 


(0) 




yi 




H, 


( 1 ) 


these equations read 


(40) 


Lf) = + Airi(y + P) . 


Using the first line to express in terms of and re-injecting the result in the second 
line yields an equation for only, namely 


r(0) 

^i+1 


An 


Y 


= Ui ( ^ + Ai ) L® + A,A,^^PLf2, 




-( 0 ) 


y2- 


for i> 1 with = 1 and = Yi -|- Ai). Using (32) and setting 


(0) 


r(0) 


y2 

w = AqAi^ P , 

we recover the well-known equation 

(41) L™=L™+»La 

for i > 1 with = L® = 1, which allows us to interpret as the generating function 
of hard pieces on a linear graph with i — 1 vertices (see Fig. |10[ ), with a weight w per piece. 
The solution of this equation is known to be (see for instancejb] Eq. (6.11)) 

1 _ y*+l 


(42) 


= 


1 


where w = -- 


1 


[l+yY l-y y + y-^ + 2 ' 

If we instead eliminate from the system (40), we obtain for the very same equation 

PW =L«+zcLW 

for i > 1, now with the initial conditions = Y and = ji = Yi (this value 

can also be read from (40) as it yields = yi(Ao -I- Ai(y + P)) = Fi(Ao -|- AiQ) = Yi). 
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We immediately deduce 

with the convention = 0. Using 
-(i)_ 1 1 


LY> = 


(43) 


(1 + y)* 1-y 
1 


(1 + y)' 


Yi(l + dy) 


we obtain 

(y(l_y* + l) + (y,_r)(l_y*)(l+y)) 
1 — ay'^^'^' 


where d = 


Yi-Y 


1-y 


and a = ^ 


y 


y2 1 + dy 


6.3. Comparison with formulas (13). Combining (39) and the explicit values (42) and 
(43), we obtain from (25) the desired expressions (17). It simply remains to show that our 
definitions for y and a of Sect. 6.2 given by (42) and (43) just above, match their definitions 
of Sect. 1^ given by (12), or equivalently (14) and ([^). If so, 0 is equivalent to ( |I^ and 
we are done. 

Using for y and a their definitions of Sect. (through ( |l4| and 0) , we obtain for 
Y = Q — P and Yi (whose value in terms of P and Q is given by ( [3^ ) the parametrizations 

{a - l)y{l - ay'^) 


Y = 


Y, = 


so that 


1 + y + ay — 6 ay^ + ay^ + a^y^ + a'^y'^ 
_ (q - l)y(l - ay^) _ 

{1 + y){l + y + ay — 6 ay"^ + ay'^ + a'^y^ + a^yY 

1 dpy 


y{l - ay) 

a = —;-^ and 

1 — ay'^ 


?/2 1 + dy 


= a 


as wanted to match the definition (43) of a of Sect. 6.2 


As for y, we use the expressions (28) of Aq and Ai to get the parametrizations 

(1 - ay^)^ 


Ar] — 


A, =- 


(1 - ay){l - ay^) 

1 + y + ay — 6 ay^ + ay^ + a^y^ + a^y^ 


(1 - ay){l - ay^) 


so that 




W = AqAi^P = — _ 

Y^ y + y +2 


as wanted to match the definition (42) of y of Sect. 6.2 


A more constructive approach consists in starting instead from the definitions of y and 
a of Sect. 6.2 (through (42) and ( |4^ ) and recovering the parametrization (14) of Sect. 
From (28) (and t, = P{1 — P — 2Q)), Aq can be expressed in terms of P and Q, as well as 
Y [= Q — P) and Yi via ([^. This leads to 


, , Y2 p(i_q_p) 

^ = ^oAi^P - - (i_2Q)2 


1 

“i 


y + y"^ + 2 


hence we deduce the (so-called characteristic) equation 
(1 - 2Q)2 - (2 + 2/ + y-')P(l - P - Q) = 0 . 
This in turn leads to 


d = 


Yi-Y 


P 


1 d + y 


P-y{l-P-2Q) 


l-P -2Q 


y^l + dy y^(l-P-2Q-yP) ' 
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Using this latter equation to express Q in terms of P, y and a, namely 
_ P{1 + y){l - ay^) 1 

2 y(l -ay) + 2 ’ 

and plugging this value in the characteristic equation above, we find that P is determined 

by 


-y{l - ay)^ + P{1 + y + ay 


from which (14) follows. 


6 ay^ + ay^ + a^y^ + a^y^) = 0 , 


7. Conserved quantities 

The explicit formulas © (resp. or equivalently ©) are typical expressions for 
the solutions of discrete integrable systems. A deeper characterization of the integrability 
of the system <§ (resp. © or @) is the existence of a number of discrete conserved 
quantities, i.e. quantities whose expression depends explicitly on some positive integer (called 
d below) but whose value turns out to be independent of this integer. In the case of bicolored 
quadrangulations, it has already been recognized that these conserved quantities may be 
obtained by looking for a direct combinatorial derivation of Fn in the slice formalism. Let us 
first recall this construction and then see how it extends to the case of our second weighting 
governed by local maxima. 


7.1. Conserved quantities for the first weighting. The slice decomposition described 
in [3 [in] applies more generally to pointed rooted quadrangulations with boundaries, that is, 
quadrangulations with a boundary, having a root-edge on the boundary oriented with the 
boundary-face on its right, and having a pointed vertex (which might not be incident to the 
boundary-face); the label d(v) of each vertex v is now the distance from v to the pointed 
vertex vq. For Q such a map, the canonical bicoloration of Q is the vertex bicoloration in black 
and white where the root-vertex (origin of the root) is black and any two adjacent vertices 
have different colors. As before, a local maximum (or “local max” for short) for the distance 
is a vertex v such that d{v) = d(v') 1 for every neighbor v' of v. For n > 1 and d > 0, let 

B'if’ be the family of admissible pointed rooted quadrangulations with a boundary of length 
2 n, where the root-vertex is at distance at most d from vq, and is one (possibly not unique) 
of the boundary-vertices that reach the smallest distance from vq. Let be the generating 
function of where each black vertex (resp. white vertex) receives weight t, (resp. to) 
except for the pointed vertex that receives weight 1. And let Zj'^( 2 n; {Bi}i>i, {Wi}i>i) be 
the generating functions of paths of length 2 n starting and ending at height d and staying 
at height at least d all along, made of elementary steps with height difference ± 1 , with each 
descending step from height i to height i — \ weighted by if i = d mod 2 and weighted 
by lUi if z = d -|- 1 mod 2 . 

Note that B^^ is nothing but the set of rooted quadrangulations with a boundary of 
length 2 n so that F^'^ = Fn. Then, as explained in 13 EQ], the slice decomposition described 


Sect. 


3.1 


for maps in Bn^ applies more generally for maps in B^'^ and yields 


Now let B^f'^ be the subfamily of B^'^ where the pointed vertex is different from the root- 
vertex, and let Fn'^'^ be the generating function for the subfamily Blf'^ where the weights are 
specified as in Fn^'^. For Q G B^'^, with vq the pointed vertex and v the root-vertex, let e be 
the first edge of the leftmost geodesic path from v to Vq. This edge cannot be a boundary- 
edge as otherwise, v would not reach the smallest distance from vq among boundary-vertices. 
We may cut along e (starting from v) so as to duplicate e into two edges Ci, e-i (with 62 before 
ei in ccw order around the new map) and duplicate v into two vertices vi,V 2 (see Fig 0 . 
Let Q be the pointed rooted quadrangulation with a boundary of length 2n -I- 4 that is 
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Figure 11. Left: a pointed rooted quadrangulation Q with a boundary. 
Right: the pointed rooted quadrangulation Q obtained by cutting along the 
first edge of the leftmost geodesic path from the root-vertex to the pointed 
vertex. 


obtained by erasing 62 , taking vi as the new root-vertex, and keeping uq as the pointed 
vertex. Denoting by d = (di,... ,^ 271 - 1 - 4 ) the distances from vq of the successive boundary- 
vertices (starting from ui) in ccw order around Q, we have the conditions that di+i = ± 1 

for i S {1, • • • ,2n + 3}, di equals the distance of v from vq in Q so that di < d, di> di for 
alH S { 1 , • • • ,2n + 1 }, d 2 n-i-i > di (indeed, by the effect of cutting along the first edge of 
the leftmost geodesic path, the distance of V 2 from vq is strictly larger than the distance of 
Vi from uq), and d 2 n +4 = di — 1. The bipartiteness of Q implies that d 2 „+i = di mod 2, so 
that the last entries of d must be (di -I- 2, di -I- 1, di, di — 1). Hence, if for A: > 1 we denote by 
{Wi}i>i) the generating function dehned as Zjj^{2n; {lFi}j>i), 

but with the restriction that the k last steps of the path are descending, then the slice 
decomposition applied to Q gives 

+ 2; {VFj,>i) • H,, 

where the factor F accounts for the (black) root-vertex being duplicated and the factor Bd 
accounts for the last descent di,di — 1. Hence for each n > 1 we have the conserved quantity 


(44) F„ = Z+^{2n- {lFj,>i) - ^Zf^'^{2n + 2; 


t. 


Bd. 


The first two conserved quantities, n G {1)2}, are (with * = d -I- 1): for all A > 1 (with 
Bo^O) 

Fi = W,- ^B,+iW,B,_u 

F 2 = ITf + H.+ilF,-^(W', + H,+i + 1F,+2)H,+i1T,H,_i. 

^9 


Shifting in (44) all path heights by —d and replacing Bi and Wi by and 14}+^ so as to 
compensate this shift, we get, upon sending d —>■ 00 the identity 

(45) F„ = Zo+o(2n; B, IF) - ^Z+;,^^i2n + 2-,B,W)-B 


t. 


(with some obvious notations) which, using the identities Z(J’g^^(2n -|- 2; B, W) = Z^Q{2n -f 
2; B, W) — Zf^Q(2n; B, W) ■ W and B = t, + B{B + 2W), is easily transformed into (27). 


7.2. Conserved quantities for the second weighting. We may now play a similar 
game for the quantities Jn to obtain conserved quantities involving the generating func¬ 
tions {Pi,Qi}i>i. Let Jn^'^ be the generating function of where each local max (resp. 
non local max) receives weight t, (resp. to) except for the pointed vertex that receives weight 
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Figure 12. Left: a pointed rooted map M with a bridgeless boundary. 
Middle: the pointed rooted map M obtained by cutting along the first edge 
of the leftmost geodesic path from the root-vertex to the pointed vertex. 
Right: the associated quadrangulation Q with a boundary (see text). 


for maps in Bn^ applies more gen- 


1. And let Z^^{2n; {Pi}i>i,{Qi}i>i) be the generating function of paths of length 2 n start¬ 
ing and ending at height d and staying at height at least d all along, made of elementary 
steps with height difference ±1, with each descending step from height i to height i — 1 
weighted by Pi if just after a descent and weighted by Q i if just after an ascent. 

Again, the slice decomposition described in Sect, 
erally for maps in and yields 

Jn ^ {Pi}i>l, {Qi}i>l)- 

Let be the family of rooted pointed general maps with a bridgeless boundary of length 

n, where the root-vertex is at distance at most d from the pointed vertex vq, and is at least 
as close from vq as any other boundary-vertex (here boundary-edges are directed ccw around 
the map while inner edges are bi-directed; the distance-label d(v) is the length of a shortest 
direc ted p ath starting from vq and ending at v). The Ambjprn-Budd bijection described in 
between Ain^ and B^'^ extends verbatim (using the same local rules, and having 


2.2 


Sect. 

the same pointed vertex and the same root-vertex in corresponding maps, see mm) to a 
bijection between B^f^ and so that is also the generating function of maps in 

with a weight t, for each non-pointed vertex and a weight to for each inner face. 

Let be the subfamily of where the pointed vertex is different from the root- 

vertex, and let be the generating function of the subfamily B^f'^ where the weights are as 
in J^\ Then the Ambjprn-Budd bijection ensures that is also the generating function 
of with a weight t, for each non-pointed vertex and a weight to for each inner face. For 
a map M G Ai)a \ let e be the first edge on the leftmost geodesic path from the root-vertex v 
to the pointed vertex vq (note that all the edges on this path are inner edges). Again we can 
cut along e (starting from v) so as to duplicate e into two edges Ci, 62 (with 62 before Ci in 
ccw order around the map) and duplicate v into two vertices ui, U 2 , and take vi as the new 
root-vertex (see Fig. 12). The map M thus obtained (as opposed to the quadrangulated case 
we do not delete 62 ) is a general map with a bridgeless boundary of length u-|-2. If we denote 
by 5i,..., Sn +2 the distances from the pointed vertex vq of the successive boundary-vertices 
(starting with ui) in ccw order around M, then equals the distance of v from vg in M so 
that (5i < d, di > for i G {1, • • • , n -|- 1 }, i5i+i < -b 1 for i G { 1 , • • • , n -b 1 }, 5n+i > (by 
the effect of cutting along the first edge of the leftmost geodesic path) and 6 n +2 = — 1. In 

particular if we reroot the map at the vertex between ei and 62 , we get a map M' G 
We may then take the image Q' G ^i+ 2 ^^ of by fbe Ambjprn-Budd bijection, and denote 
by Q the quadrangulation with boundary obtained from Q' by shifting the root position 
by one in ccw order around Q'] note also that the number of local max (resp. non-local 










24 


ERIC FUSY AND EMMANUEL QUITTER 


max) of Q equals the number of inner faces (resp. the number of vertices) of M', which 
is also the number of inner faces (resp, the number of vertices plus 1) of M. Let again 
d {dij • ■ -, 1 ^ 271 + 4 ) be the distances from the pointed vertex vq of the successive boundary- 
vertices (starting with rii) in ccw order around Q. By the local rules of the Ambj0rn-Budd 
bijection, d is obtained from the sequence ^i, ..., (5„_|_2 where for each z S {1, • • • , n -b 1}, we 
insert between Si and (5i+i the subsequence (of length Si — 6i+i -1-1) Si + l,di,..., Si+i -I- 1. 
It is then easy to check that d satisfies the following conditions: di = Si, di+i = diztl and 
di > i5i for z G {1, • • • , 2n-|-3}, d 2 n +4 = — 1, and d ends with di -b 2, di -b 1, di, di — 1 (since 

Sn+i > di). Hence, if for fc > 1 we denote by Z'^'^"^{‘2.n-,{Pi\i>i,{Qi\i>i) the generating 
function defined as Z^^{2n; {Pi}i>i, {(5i}i>i), but with the restriction that the k last steps 
of the path are descending, then the slice decomposition applied to Q gives 

+ 2; {Pj.>i, {QJ,>i) • Pd, 

where the factor A accounts for the root-vertex of M being duplicated, and the factor Pd 
accounts for the last descent di, — 1. Hence for each rz > 1 we have the conserved quantity 


(46) 


Jn — ^ddi'^'^’ {Qi}i>l) ~ {Pi\i>l^ {Qi}i>l) ' 


t. 


Pd- 


Remarkably this has exactly the same form as the bicolored conserved quantities (|44| , up to 
changing {Pi,Qi}i>i for {Bi,Wi}i>i and taking the “hat” variants of the path generating 
functions. The first two invariants, n G {1)2}, are (with z = d — 1): for all z > 1 (with 
Ho = 0) 

Jl — Qi ^Qi+lPiPi—1, 

J 2 = QI + Qi+lPi — —[[Qi + Qi+l)Qi+l + Qi+2Pi+l)PiPi-l- 

^9 


As before, upon sending d —> cx) in (461, we get the expression 
Jr, = 4%(2zz; P, Q) - ^Zo+o"^(2n + 2; P, Q) • P 

^9 

(with straightforward notations). Upon using P — B, Q = W and comparing with (45), this 
provides another (computational) proof of the identity J„ = Fn by noting that Z()"Q(2n; P, Q) = 
Z(}o(2n;P, Q) and Z(^Q^^(2n -b 2;P,Q) = ZQ_Q^{2n + 2;P, QjH Finally, from ([^, we get 
Hi = (Qi —Pi) = to—t. + (Qi—Pi)Qi, hence Yi = (to—t,)/(l—Qi). Using the first conserved 
quantity above, we deduce Qi = Ji = Q — QP^/t, so that Yi = t,(to — t,)/(t, — t.Q-bQP^) 
which upon expressing t, and to in terms of P and Q via (15), reproduces the expression 
(321) for Yi. 


8. Conclusion 

In this paper, we presented a comparative study of two statistical ensembles of quad- 
rangulations. We first showed how the corresponding slice generating functions {Bi,Wi for 
the first ensemble and Pi,Qi for the second) appear as coefficients of the same quantity 
F{z) = J{z), expanded as a continued fraction in two different ways. The slice generating 
functions may then be written as bi-ratios of Hankel-type determinants and explicit formulas 

^ The identity ZgQ{2n;P,Q) = Z^Q(2n;P,Q) is easily proved by noting that the equation Z = 1/(1 — 
z(Q — P) — zP Z which determines the generating function Z = 'JZn>o Q)z" is identical to that, 

Z = 1/(1 — zQ/(l — zP Z)) which determines the generating function Z = Z^Q(2n; P, ( 5 ) 2 :^, hence 

Z = Z. The identity Z^q'^ (2n.+2; P, Q) = Zqq'^ (2n.+2; P, Q) follows by noting that Zqq'^ (2n+2; P, Q) = 
Z^Q(2n+2; P, Q)—Z^Q(2n; P, Q)-Q and similarly (2n+2; P, Q) = Z^Q(2n+2; P, Q)-Z^Q(2n; P, Q)-Q. 
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Figure 13. A example of heap of 7 pieces sitting on top of the graph 
Q of Fig. 14 with base {1,8} (we indicated in light bine the “shadow” of 
those pieces which can move freely and hit the vertices of the graph). The 
diameter of the pieces is adjnsted so that pieces sitting on top of vertices 
which are adjacent in Q cannot pass through each other. 


may be obtained, at the price of some conjectnred expression for some intermediate quantity 
in the second ensemble. 

To conclude, we would like to emphasize that our two ensembles may be viewed, in some 
sense, as the two extremal elements of a very general family of statistical ensembles as 
follows: by definition, the second ensemble gives a particular weight to those vertices which 
are local maxima for the distance to the root-vertex. Similarly, the first ensemble may be 
viewed as the ensemble which gives a particular weight to those vertices which are local 
maxima for the distance to the root-vertex modulo 2. Indeed, this distance modulo 2 is 0 for 
black vertices (recall that the root-vertex is black) and 1 for white vertices so that all white 
vertices are local maxima. In this respect, note also that performing the passage from the 
quadrangulation to the general map in the bijection of Fig. may be viewed as applying the 
Ambj0rn-Budd rules, taking as labeling the distance modulo 2. 

Denoting by d{v) the distance from a vertex v to the root-vertex in a rooted quadrangu¬ 
lation with a boundary, we may more generally consider statistical ensembles which give a 
particular weight to those vertices which are local maxima for some labeling i(d(y)), with 
d !->■ l{d) some given function. Without loss of generality, we may set .^(0) = 0 and, if we 
wish to apply the Ambj0rn-Budd rules to transform our quadrangulation into a general map, 
we need that \i{d) — £{d — 1)| = 1 (it also seems natural to impose that £{d) remains non 
negative so that the root-vertex cannot be a local maximum). It is likely that slice generating 
functions in this ensembles may appear as coefficient of F(z) = J{z), once expanded as a 
continued fraction with some appropriate structure, being a mixture of the Stieljes-type and 
of our new encountered type. At this stage, it is interesting to notice that, in their study 
of finite continued fractions HIE], Di Francesco and Kedem introduced precisely a whole 
family of such “mixed” fractions as well as some passage rules on their coefficients to go from 
one to the other without changing the actual value of the fraction. It is very tempting to 
speculate that their study may be extended to infinite continued fractions to describe our 
more general ensembles. 
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Figure 14. The semi-infinite graph Q and the finite graph 


Appendix A. A proof of the formulas (251 and (261 


As in our proof of formulas ( [2^ and (|2 m is based on the theory of heaps of pieces. 

The reader is invited to consult m for the basics of this theory. 

Let us simply recall what we mean by a heap of pieces on a graph Q, supposedly connected, 
planar, and drawn in a horizontal plane for simplicity. Imagine to complete the graph by a 
set of vertical half-lines, with a half-line starting from each vertex of the graph. Informally 
speaking, a heap is a collection of pieces threaded along these half-lines. Each piece therefore 
sits on top of a given vertex and may move freely along the corresponding vertical half-line 
as long as it does not meet another piece. More precisely, the pieces are supposed to be 
designed so that two pieces may not pass each other if they sit on top of the same vertex or 
if they sit on top of adjacent vertices. 

Given a subset B of the set of vertices of Q, a heap of pieces is said to be of base B if, 
moving its pieces as far as possible to the bottom of the half-lines, the set of those vertices 


hit by a piece forms a subset of B (see Fig. 131. 

A fundamental remark is that, from the relation ([^, J{z) may be viewed as the generating 
function for heaps of pieces on the semi-infinite graph Q of Fig. |14[ with a weight zYi per 
piece sitting at position i along the graph, and whose base is {1,2}. Similarly, we may 
interpret J{z) as the generating function for the very same heaps, but now with a weight 
zYi per piece sitting at position i. Let us finally introduce the quantity 

K{z) = 1 -I- z Yi J{z) = - 


1 - z- 


Y 


1-z- 


Yo 


l-zY.- z- 


Y4 


1-zY,- 




1 - 


which is the generating function for heaps of pieces on the graph Q again with a weight zF} 
per piece sitting at position i along the graph, but now with base {!}. 

From the definition (23) of the ’s, we have 




n>0 


J{z) = 


n>0 


J-nZ 


so that all the j„’s have a direct interpretation as enumerating heap configurations made of 
|n| pieces. 
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Let US now consider the analogs K^°‘\z) and J^°‘\z) of J{z), K{z) and J{z) 

respectively, viewed as heaps generating functions, now defined on the finite graph of 
Fig. In other words, we set 




1 


1 - zYi - z- 

I-ZY3- z. 


1^2 


Yi 


1 - zY2a-3 - 


Yon-^ 


1 — Z I 2 Q-I 




1 


1 - zFi - z- 


1^2 


1-zY.- z- 


>4 


i^(“)(z) = l + zYi j(“)(z) 


1 — ^^2a-3 ~ 


^20-2 


i — z I2Q-1 

1 


1 - z- 

1 - z 



1 - zFa - 


Yi 


We finally define the analogs of j„ via 






J(“)(z) = ^£ 




n >0 


n >0 


1 - zy2a-3 - ^ 


F2a-2 
'1 - zF2a-l 


so that (n > 0) enumerates heap configurations of n pieces on with weights Yi and 
base {!}, and (n > 0) enumerates heap configurations of n pieces on with weights 
Yi and base {1,2}. 

It is now a standard result of the theory of heaps of pieces [12] tha10 


(47) 


^(a). ^ ^ X^‘-\^,-zY2,-zY3,--- ,-zy2o-l) 

x(“)(-zyi,-zi2,-^F3,-- - ,-zr2c-l) 
J(a), . ^ X(")(0,0,-zf3,--- ,-^^ 20 - 1 ) 

Xi<Y{-z%,-Z%,-zY3,--- ,-zy2a-l) 


^ A sketch of the proof is as follows: given 3 , consider pairs made of a heap configuration "H of 

base B together with a configuration 'HP of hard pieces, drawn on top of the heap. For such a pair, let S be 
the set of pieces that can be moved up freely to infinity, and when pushed downward either are blocked by 
a piece (that has to be in H) or hit a vertex of the base B. Consider the following transformation: if S is 
not empty, pick the piece p £ S oi smallest index and change its status (from H to HP if p £ H, from HP 
to H if p £ HP); if S is empty do nothing. This transformation is easily seen to be an involution (which 
leaves S invariant), and, if we assign a weight 2 per piece in the heap and —2 per piece in the configuration 
of hard pieces, the weight is multiplied by —1 for each configuration which changes under the involution. 
The generating function for the pairs, which is the product of the generating function for heaps with a 
weight 2 per piece times that of configuration of hard pieces with a weight —2 per piece, therefore reduces 
to those pairs for which S is empty. It is easily seen that this situation corresponds to an empty heap and a 
configuration of hard pieces made of pieces which do not belong to B. The corresponding generating function 
is nothing but that of configurations of hard pieces with a weight —2 per piece not in B and 0 per piece in B. 
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Figure 15. The maximally occupied configuration of hard pieces on the 
graph made of a pieces (here represented by gray circles) sitting on 

all the sites with an odd index. 

where X^°\yi,y 2 ,y 3 , ■ ■ ■ , y 2 a-i) denotes the generating function of hard pieces on the graph 
Q^°‘\ each piece sitting at position i receiving the weight Recall that, by definition, in a 
configuration of hard pieces, each vertex of the graph is occupied by at most one piece, with 
no two adjacent vertices occupied simultaneously. Note that the positions of the O’s in the 
numerators correspond to the location of the vertices of the corresponding base of the heaps 
({1} and {1, 2} respectively). Clearly, on the graph we can put at most a hard pieces. 
Moreover, this maximal situation is achieved by a single configuration with all sites with 
odd index occupied (see Fig. @). The quantity zYi, —ZY 2 , —zY^, • • • , —zY 2 a-i) is 

therefore a polynomial of degree a in z that we write 

Ct 

X(‘^'>i-zY,,-zY2,-zY3,- ■ ■ ,-zY 2^.,) = J2i-^r^tHyuY2,Y3,- ■ ■ ,Y2a.-l) 

m—0 

where Xm \yi,y 2 ,y 3 , • • • j 2 / 2 a-i) denotes the generating function of exactly m hard pieces 
on the graph each piece sitting at position i receiving the weight yi. Clearly, both 

X(‘^)(0,-zY2,-zY3,--- ,-zT 2 „_i) and X(“)(0,0,• • • ,-zYj^.i) are polynomials of 
degree a — 1 in z. 

Let us now come to our fundamental identities. We have 

v(a) _ v(a) T>(a) 

( 0 ) = ( 0 , 0 )) 

with the short-hand notations 

= Xi“)(Fl, I 2 , is, • • • , i^2a-l) 

Xi“)(0)=XM(0,y2,T3,-- - ,l2a-l) 

X^'> =X^^\YuY2,Y3,--- ,F2o-i) 

li“)(o,o) = xi“)(o,o,y3,--- ,y2a-i) • 

To explain these identities, let us analyze the structure of a conhguration C of hard pieces 
on In C, a number k of pieces occupy even sites 2ji, 2j2, • • • , 2jfe with 1 < < j 2 < 

••■< jk < ct — 4 and — ji > 1 for £ = 1, • • • ,k — 1. The set of available odd sites is 

Odd = {1, 3, 5, • • • , 2a-l}\{2ji- 1, 2j2-l, 2j2 + l, • • • , 2jfe-l, 2jfc-|-l} and satisfies 

|C>dc?| = a — 2k. A number k' of pieces occupy a subset {2j[ — 1, 2j2 — 1, • • • , 2j(,, — 1} of this 
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Figure 16. Left: an example of configuration C of hard pieces (represented 
in gray) on . Right: the associated configuration C of hard pieces 
(represented in light blue), obtained by keeping the particles sitting on even 
vertices and, in the ensemble of odd vertices which are not adjacent to the 
occupied even vertices, exchanging the occupied and un-occupied sites. 


set. In Xm^ (corresponding to a situation where k + k' = m), any occupied site i receives the 
weight Yi, so the weight of the conhguration C is Y 2 j^Y 2 j 2 ■ ■ ■ Y 2 j^ x F' 2 jj_iy 2 j'-i • ■ ■ 

Let us now consider instead the configuration C where, again, the sites 2ji, 2j2, • • • , 2jfc are 
occupied but now the complementary of {2j[ — l,2j2 — 1, - • • , 2j(., — 1} in Odd (namely 
Odd \ {2j( — 1, 2j2 — 1, • • • , 2j(., — 1}) is covered by pieces. Clearly, going in from C to C 
provides a bijection between configurations C with k -\- k' = m pieces and configurations C 
with fc + (a — 2k) — k' = a — m pieces. In Xa-m, the configuration C receives the weight 


Y2nY2,. 


■Y, 


n 

iGOdd 


‘^jk 


Y2nY2n 


■Yn 


2jk 




■Y, 


2jL-1 


-1 ^ 2 + 1 ^ 2 ^ 2-1 > 2^2 +1 


■ Y2j^-iY2j^ + i 


"" n 

i^Odd 

Y2j^Y2j^ ■ ■ ■ Y2j^ X Y2f_^-iY2f^-i ■ ■ ■ Y2 j' I 


Y2nY2j, 


_ 

Y1Y3Y5 ■ ■ ■ Y2o,-1 

Y2j^ X Y2f^-iY2f^- 


■Y2jI,-1 


X, 




since = YiYaYs • • • F 2 q,_i. From the bijection C i-A C, we therefore deduce immediately 
the first equality in (481. To get the second equality, we note that Xm\o) enumerates 
conhgurations C with m pieces such that the site 1 is not occupied by a piece. Two situations 
may then occur: either site 2 is occupied or not. In the first case, the bijection C ^ C will 
generate a configuration C where site 2 is occupied (and site 1 does not belong to Odd) while 
in the second case, it will generate a configuration where site 2 is empty and site 1 (which 
belongs to Odd) is necessarily occupied (since it was empty in C and the empty and occupied 
sites get exchanged in the bijection for those odd sites belonging to Odd). To summarize, 
in the configuration C, either site 1 or site 2 must be occupied. The restriction of xj^2m to 
these configurations yields X^-m ~ 0)) hence the second equality. 
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From (48), we deduce 


Id ^2 Fa >da-l 


y(«) ••• - 

I 5 5 5 5 

Z Z Z Z 


dda-l 


x(“) 

\ z z 


-x(“)(0,0,-zy3,... ,-zf2„_i)) 


and therefore, taking the ratio of the two lines and using (47), 

1 + y Q = 1 - j(“)(z) . 

The finite continued fraction case of Sect. 5.3 corresponds precisely to a situation where 
J(z) = J^“^(z) and J{z) = J^°‘\z). The above formula explains the first identity in (26) 
while the second identity is guaranteed by the relation J^°‘\z) —)■ 1 when z —)■ 0. This 
concludes the proof of (261. 

We now prove (25) by computing explicitly the determinants Hf*'^ and in terms of 
the id’s. More precisely, let us show that, for i > 1, 


(49) 





jo 

Hf> ^ 



jl 


jo 

jl ■ 

■ j^- 

-1 


j-i^-2) 


■■ jl 





J2 

= 


jl 

J2 • 

■ j^ 



»-2 


^21-4 

h2i-3 


Y^ 

hdi-l 


= w Fa Fs • • • F2 ,._i H. 


2 i-l 


(0) 


Once these formulas are proved, the relations (25) indeed follow immediately. 

A first crucial point is the existence of a linear relation between the ’s, namely 


(50) ^ (-l)™Ai“) = 0 for all integers n . 

m=0 


Indeed, writing the first identity in (47) as K^°‘\z)X^°‘'>{—zYi,—zY 2 , —zFa, • • • , —zF 2 q_i) = 
X(“)(0, —zF 2 , —zY^, • • • , —zF 2 a_i) and extracting the term of order z", we immediately see 
that ( [^ holds for any positive integer n> a since A^“)(0, —ZF 2 , —zFa, • • • ,—zF 2 q_i) is 
a polynomial of degree a — 1. Similarly, writing the second identity in (47) as J^°‘\z) x 
X(“)(—zFi, —ZF 2 , • • • ,—zY 2 a-i) = A(“)(0,0, —zla, • • • ,—zh^a-i), a polynomial of degree 
a — 1, we find that 


a 

E 

m —0 


I — ! y(“) _ n _ ( \ ' /I y(“) •(“) 


Ad 


m'—O 


for n' > a. Here we have set m' = a — m and used Setting n = 


a — n' < 0, we deduce that (50) also holds for any non-positive integer. It remains to show 
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that it is valid in the range 1 < n < a — 1. For n in this range, we have 

an a 

E I — I i™ Y(“) —1 ("—I 1"* V(“) o(a) 

\ Jn—m / j \ Jn—m ' / j \ Jn—m 

m—0 7n—n-\-l 

a—n— 1 

= (-l)"X“(0)+ (-!)“ ^ 


m—0 


p—0 


= (-irx:(o) + (-ir4“) ^ . 

p—O 


Now since l<a — n<a — 1, we also have 

a—n 


p—0 

so that we eventually get 


O' 

= (-1)" (x“(0) - - li“l(0,0))) = 0 . 


m—0 


The linear relation (50) therefore holds for all integers n, as stated. 

Let us now come to the computation of and Since jn, n > 1, enumerates heaps 
of n pieces on the graph Q with base {!}, the pieces cannot reach sites with index more than 
1 for n = 1 and 2n — 2 for n > 2. In other words, jn enumerate heaps which “live” on 
therefore on for all n < i — 1. As for jn, n < —1, it enumerates heaps of \n\ pieces 

on the graph Q with base {1,2} so that the pieces cannot reach sites with index more than 
2|n|, therefore “live” on therefore on for all |n| < j — 2. In other words, we 

have 

jn = j^~^^ for n = 0,1, 2, • • • ,i-l 
j-n = j-n^^ for n = 0,1, 2, • • • ,i-2 

In the determinant the only term which does not “live” on is and it is 

easily seen that 

with an additional term corresponding to the unique heap that hits position 2i — 2. Using 
the linear relation (50) for a = i — 1, we may thus rewrite as 


i7io) = 


■ ■ ^21-2 

j-ii-2) ■ ■ 

jo 

0 

j-{i-3) ■ ' 

jl 

0 

jl 

■ jr- 


= ■ ■ ■ F 2._2 . 


Similarly, in the determinant the only term which does not “live” on is ji and 

it is easily seen that 


J9=jf +Yi{Y2Yi---Y2^-2) 
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with again an additional term corresponding to the unique heap that hits position 2i — 2. 
We may thus rewrite as 




jo 

jl 


Jo 


J,_l Yi{Y2Y4---Y2^-2) 


= Fl (1-2^4 • • • r2*-2) K 


(0) 


Combining the two above formulas and replacing the l^’s by their value in terms of the l^’s, 
we deduce the recursion relation 

2 


.Y 3 Y 5 


Y 2 Y 4 


Y2i-4 

Y2i-3 


^2z-2 ^(0) 
^22-1 


i-2 


(0) _ 1 w(0) = = (^ 2 /^ 3 ). 


for i > 3 with initial conditions = 1 and H 2 
The first line of eq (49) follows immediately. As for the second line, it follows from 

r(0) 




H. 


(0) 


= yi(y2r4---r2,_2)- 


H, 


(0) 


= ri(r2i4' 


•1^2.-2)(^^ 


^21-1 

b2i-2 


= Y^Y^Yz ■ ■ ■ ^ 2.-1 


The ab ove d erivation of Eq. (|25|) extends verbatim to the case of the finite continued fraction 
of Sect. 


5.3 


by limiting to i < a — 1 the range of allowed values for the index i in and 
This range is precisely what is needed to compute 11,12, • • • , lla-i- 

Appendix B. A proof of the formulas (|^-([M|) 

The quantities Z{z\P,Q), Z{z;P,Q) are specializations of J{z) and J{z) (viewed as de¬ 
fined from the Yjs through their continued fraction expansions) to the case where 


l2z-l = E , 


Y 2 , = P 


12.-1 = 


y2^ = ^ = p 


for all i > 1. Consequently, Zn, Zn and kn are the corresponding specializations of J^, 
Jn and jn- The analysis of Appendix A applies to arbitrary ll’s. In particular, kn has a 
direct interpretation in term of heaps of pieces on the graph Q for all n. For n > 0, kn 
enumerates heaps of n pieces of base {!}, with weights Y and P for pieces on odd or even 
sites respectively. For n < 0, kn enumerates heaps of |n| pieces of base {1, 2}, with weights 
Y and P for pieces on odd or even sites respectively. Let us thus introduce the quantities 
(analogs of corresponding to a restriction of the heaps to the graph From 

(50), we deduce immediately 

Oc 

(— kl^2ni = 0 for all integers n . 

m—0 

where = xin\y^ P, F, P, • • ■ , F) is the generating function of configurations of exactly 
m hard pieces on the graph Setting a = i — 1, this equation reads equivalently 


m—0 




i-iy 


.X', 


i-l) 

-1-77, 

d»-l) 


= 0 for all integers n . 


Now, from their heap interpretation, it is clear that kn = k. 


(7-1) 


for n = 0,1, • • • , i — 1 


as well as for n = —1, —2, • • • , —{i — 2). For 2 < n < i, all the kp ^^’s appearing in the 
above formula may thus be replaced by fcp’s and (35) follows. For n = 1, the only term 
which gets out of the graph is for m = 0, since (the two 

indeed differ by the contribution of the heap made of one piece on each even site from 2 to 
2(j — 1)). This explains the right hand side = p*-i/p2(7-i) (36). 

For n = i + 1, the only term which gets out of the graph is for m = z — 1, since 

















COMPARING ENSEMBLES OF QUADRANGULATIONS VIA CONTINUED FRACTIONS 


33 


2 a—1 2a 


3 

1 




2 a-2 


4 

2 


Figure 17. The graph 


= ki — YP^ ^ (the two indeed differ by the contribution of the heap made of one piece 
on site 1 and one piece on each even site from 2 to 2(i — 1)). This explains the right hand 
side = (—in the second line of (36). 

Eqs. (37) and (38) can be proved in the same way but their proof now relies on a restriction 
of the heaps to the graph of Fig. 


17 


_ Let us first analyze the heap generating functions 

on this graph: denoting, for n > 0, the generating function of heaps of n pieces on 
of base {1} and with weight Yi per piece sitting on site i. It is easily seen that 
also corresponds to enumerating heaps of n pieces on the graph provided we assign, 
instead of Y 2 a-i, a weight Y^^-i = Y 2 a-i + Y 2 a to pieces sitting at position 2 q: — 1. The 
same remark holds for configurations of m hard pieces on which are enumerated by 
Xm\Yi,Y 2 , ■ ■ ■ ,¥ 20 - 2 ,T 2 q,_i). In order to use directly our previous results (obtained for 
we are thus lead to define for n < 0 as enumerating heaps of n pieces with base 
{1,2} on the graph with weights Yi built via the same expression (21) as before with 
F 2 a-i replaced by Y^^-i- Getting back to for n < 0 therefore enumerates heaps 

of n pieces with base {1,2} on this graph, with weights Yi as in (21) for pieces on sites 
1,2,3, • • • , 2 q! — 3 and the special weights ^ 20 - 2 / (Y 2 a- 3 (Y 2 a-i + Y 2 a)) for pieces on the site 
2 q: — 2, l/(F 2 a-i + Y 2 a) for pieces on the site 2a — 1 and 0 for pieces on the site 2a. With 
this definition, we have the analog of (50), namely 


(-1)™X((“) = 0 for all integers n . 

m—0 


where enumerates configurations of m hard pieces on ^d“)_ 

Let us now specialize this result upon introducing, for n > 0, the generating function kn°‘'^ 
of heaps of n pieces on ^ of base {1} and with weights Y and P for pieces on odd or even 
sites respectively. From the above discussion, for n < 0, must be defined as enumerating 
heaps of n pieces on with weight Y for pieces on odd sites 1,3, 5, • • • , 2a — 3, P for 
pieces on even sites 2,4, • • • , 2a — 4 and the special weights P/{Y(Y + P)) for pieces on the 
site 2a — 2, l/iY + P) for pieces on the site 2a — 1 and 0 for pieces on the site 2a. With 
these definitions (and A;q“^ = 1), we have 

Ol 

(— 1 )™ 2 ;(G) — 0 for all integers n 

m=0 
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where x'm^ enumerates configurations of m hard pieces on with weight Y (resp. P) 
per piece sitting on an odd (resp. even) site (in particular Xq°‘^ = 1). Setting a = i — 1, the 
above equation becomes 


i-l 

E 

m—O 


(-ir 


^n-m^ — 0 for all integers n. 


Again, from their heap interpretation, it is clear that = k- 


'b-i) 


for n = 0,1, 


, i as well 


as for n = —1,—2, •• • ,—(i — 3). For 2 < n < z, all the fcp* ^^’s appearing in the above 


formula may thus be replaced by fcp’s and (37) follows. For n = 1, the only term for which 
this substitution fails is for to = z — 1 , since — — P/ (f^(^ + -P))) (the 

two indeed differ by the contribution of the last piece, at position 2(z — 2), in the the heap 
made of one piece on each even site from 2 to 2(z — 2)). This explains the right hand side 


(_!)* ip* ^(P — Pf(Y{Y + P)))x'^\^'^ = (—!)* ^ in the first line of (38) (note 


that 


{Y + P)). For n = z + 1, the only term for which this substitution fails is 

for TO = 0, since = h+i — YP^~^{Y + P) (the two indeed differ by the contribution 

of the two heaps made of one piece on site 1, one piece on each even site from 2 to 2(z — 1) 
and a last piece at position 2z — 1 or 2z). This explains the right hand side YP^~^{Y + P) 
in the second line of (38). 
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